METASTABILITY FOR SCALAR CONSERVATION 
LAWS IN A BOUNDED DOMAIN 
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Abstract. The initial-boundary- value problem for a viscous scalar conservation 
law in a bounded interval / = {—I, i) is considered, with emphasis on metastable 
dynamics, whereby the time-dependent solution approaches its steady state in an 
asymptotically exponentially long time interval as the viscosity coefficient e > 
goes to zero. A rigorous analysis is used to analyze such slow motion for solutions 
belonging to a cylindrical neighborhood of a one-parameter family of approximate 
solutions {[^^(s 0}{£/ J to be considered as an approximate invariant manifold for 
the problem. The description consists in an ODE for the parameter ^, describing 
the position of an internal transition layer, coupled with a PDE describing the 
evolution of the distance vector v :— u — {■;£,) from the manifold. By use of 
the properties of the linearized operator at , we estimate the size of both layer 
location ^ and distance vector v. 

Key words. Metastability, slow motion, spectral analysis, viscous conservation 
laws. 

1. Introduction 

Metastability is a broad term describing the existence of a very sensitive equi- 
librium, possessing a weak form of stability/instability. Usually, such behavior is 
related to the presence of a small first eigenvalue for the linearized operator at the 
given equilibrium state, revealed at dynamical level by the appearance of slowly 
moving structures. Such circumstance comes into view in the analysis of differ- 
ent classes of evolutive PDEs, and it has been object of a wide amount of studies, 
covering many different areas. Among others, we emphasize the explorations on 
the AUen-Cahn equation, started in O HO], and the investigations on the Cahn- 
Hilliard equation, with the fundamental contributions [25l[T]. The analysis has been 
continued by many other scholars by means of a broad spectrum of techniques, 
and extended to a number of different models such as the Gierer-Meinhardt and 
Gray-Scott systems (see [29]), Keller-Segel chemotaxis system (see [HI [26]), general 
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gradient flows (see [21]) and many others. The number of references is so vast that 
it would be impossible to mention all the contributions given in the area. 

A pionereeing article in the analysis of slow dynamics for evolutive parabolic 
partial differential equations has been authored by G. Kreiss and H.-O. Kreiss [H] 
and concerns with the scalar conservation law 

(1.1) dtu + dj{u)=edlu, 

with the space variable x belonging to a one-dimensional interval / = (— £, i), i > 0. 
The primary prototype for the flux function / is given by the classical quadratic 



formula f{s) = \ s , so that equation (1.1) becomes the so-called (viscous) Burg- 



ers equation. The parameter e > is considered as a strictly positive and small 



parameter. Equation (1.1) is complemented with Dirichlet boundary conditions 

(1.2) u{-i,t)=u. and u{i,t) = u+ 

for given data to be discussed in details. 

Burgers equation is considered as a (simplified) archetype of more complicate 
systems of partial differential equations arising in different fields of applied mathe- 
matics. Inspired by the equations of fluid-dynamics, the parameter e is interpreted 
as a viscosity coefficient and the main problem is to identify and quantify its role in 
the emergence and/or disappearance of structure. 



Formally, in the limit £ — )■ O"*", the parabolic equation (1.1) reduces to a first-order 
quasi-linear equation of hyperbolic type 

(1.3) dtu + dJiu) = 0, 

whose standard setting is given by the entropy formulation, hence possessing possibly 
discontinuous solutions with speed of propagation s given by the Rankine-Hugoniot 
relation 

(1.4) sM = lf{u)] 

(where |.] denotes the jump) together with appropriate entropy conditions. In ad- 



dition, the treatment of the boundary conditions (1.2) is much more delicate with 
respect to the parabolic case, because of the eventual appearance of boundary layers. 
Focusing on the more significant example, let us assume 

(1.5) f"is) > Co > 0, /'(n+) < < /'(«_), /(m+) = /(n_), 

for some positive constant cq. The last two assumptions guarantee that the jump 
from the value ti_ to m+ is admissible and its speed of propagation, dictated by ( |1.4[ ), 



is zero. In this case, equation ( |1.3[ ) possesses a one-parameter family {f^i,yp('! 0} 



stationary solutions satisfying the boundary conditions ( 1.2[ ), given by 



where Xi denotes the characteristic function of the set I. The dynamics determined 



by initial-value problem for (1.3)-(1.2) is very simple: it is possible to prove that 
every entropy solution converges in finite time to an element of the family {U^_^^^{-; ^)} 



METASTABILITY FOR SCALAR CONSERVATION LAWS 



3 



(see Section |5j). Hence, at the level e = 0, there are infinitely many stationary 
solutions, generating a "finite-time" attracting manifold for the dynamics. Note 
that, at the hyperbolic level, there is no way of distinguishing one element from any 
other in the family of steady states {U^^^^{■] C,)} in term of stability properties. 

For e > 0, the situation is different. Apart from the well-known smoothing effect. 



the presence of the Laplace operator at the right hand-side of (1.1), together with 



the boundary conditions (1.2), has the effect of a drastic reduction of the number 
of stationary solutions: from infinitely many to a single stationary state (see [Ti]). 
Such solution, denoted here by [7^ = If^ (x), converges in the limit e — >■ 0"*" to a 

' par par v / ' ^ 

specific element of the family {U^^^^{-; ^)}. As an example, in the case of the Burgers 
equation, /(s) = | s^, there holds 

f/^ (x) = —K tanh f^— 
\2e 



par ' 



where k = Kc{e/i,U-) is implicitly defined by the relation 

K{e, u) tanh {K{e, u) /2 9) = u. 

In the limit e — > 0"*", the single steady state U^^^ converges pointwise to f/^^p := 
f/^yp(-;0). Therefore, the element of the one-parameter family {U^_^^^{■^, C,)} corre- 
sponding to ^ = exhibits a form of structural stability which is not shared with 
any other element of the same family. 

A deeper understanding of the problem can be gained by analyzing the dynamical 



properties of (1.1) for initial data close to the equilibrium configuration by means 



par 



of the linearized equation at the state U'^ 
^ p 

dtu = CsU := ed^u + dx{a{x)u) with a{x) := —f\U^^^{x)). 

In [H] it shown that, in the case of Burgers fiux /(s) = \ s^, the eigenvalues of C^, 
considered with homogeneous Dirichlet boundary conditions, are real and negative. 
Moreover, for f{u+) = f{u^), there holds as e — ?■ 

Xl = 0{e~^/') and < < VA; > 2 

for some cq > independent on e. Negativity of the eigenvalues implies that the 
steady state U^^^ is asymptotically stable with exponential rate; while the precise 
distribution shows that the large time behavior is described by term of the order 
e'^'i* and, as a consequence, the convergence is very slow, when e is small. 

Such is the picture relative to the behavior determined by an initial data close 
to the equilibrium solution tJ^ . The next question concerns with the dynamics 

■'■ par ■'■ 

generated by initial data still presenting a sharp transition from u~ to localized 
far from the position of the steady state U^^^ . Figure ll represents the solution to 
the initial-value problem for (1.1) with boundary conditions ( 1.2[ ). Starting with 



a decreasing initial datum, a shock layer is formed in an (9(l)-timescale, so that 
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Figure 1. The solution to viscous Burgers equation (1.1) with e = 0.04, 
u± = =f1 and initial datum u(x,0) = uo{x) := l/2x^ — x — 1/2. 

the solution is approximately given by a translation of the (unique) stationary so- 
lution of the problem. Once such a layer is formed, it moves towards the location 
corresponding to the equilibrium solution and this motion is exponentially slow. 

Heuristically, since, for e = 0, there exist infinitely many steady states with a 
single jump located at an arbitrary point in the interval {—£,£), for e > the 
corresponding profiles -even not being stationary solution- still solves approximately 
the evolution equation. As a consequence, two different time scales emerge: for short 
times, the solution becomes closer to a monotone transition connecting the boundary 
data, close to a space-translation of the single steady state; for long time, the profile 
slowly moves toward the equilibrium configuration. 

To describe the dynamics generated by such an initial configuration and to de- 
termine a detailed description of the relation between the unviscous and the low- 
viscosity behavior, it is rational: 

- to build up a one-parameter family of functions {f/^(-;0} such that U^{-]^) = 
U^^^ for some ^, and with the additional property that [/'^(■; ^) — )■ U^_^^^{-;^) as e — )■ 0, 
in an appropriate sense; 

- to describe the dynamics of the viscous scalar conservation law in a tubular 
neighborhood of the family {?7^(-;^)}. 

Among others, such problem has been examined in [27] and in [16] , where different 
approaches have been considered. The former is based either on projection method 
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or on WKB expansions] the latter stands on an adapted version of the method of 
matched asymptotics expansion. The common aim is to determine an expression 
and/or an equation for the parameter ^, considered as a function of time, describing 
the movement of the transition from a generic point of the interval {—£,£) toward 
the equilibrium location ^. In both the contributions, the analysis is conducted at 
a formal level and validated numerically by means of comparision with significant 
computations. A rigorous analysis has been performed in [7] (and generalized to 
the case of non convex flux in jSj), where one-parameter family of reference functions 
is chosen as a family of traveling wave solution to the viscous equation satisfying 
the boundary conditions and with non-zero (but small) velocity. The approach is 
based on the use of such traveling wave to obtain upper and lower estimates by the 
maximum principle, from which rigorous asymptotic formulae for the slow velocity 
are obtained. 

Slow motion for the viscous Burgers equation in unbounded domains has been also 
considered in literature. In [28j, it is analyzed the case of the half-line (0, +oo) for 
the space variable x, with constant initial and boundary data chosen so that speed 
of the shock generated at a; = would be stationary for the hyperbolic equation. 
The presence of the viscosity generates a motion of the transition layer, which is 
precisely identified by means of the Lambert's W function. Later, the (slow) motion 
of a shock wave, with zero hyperbolic speed, for the Burgers equation in the quarter 
plane has been considered in [12], where it is shown that the location of the wave 
front is of order ln(l + 1); the same result has been generalized in |23j in the case of 
general fluxes (for other contributions to the same problem, we refer also to [171 130] ) . 

The case of the whole real line has been examined in [13] with emphasis on the 
generation of A^— wave like structures and their evolution towards nonlinear diffu- 
sion waves. The analysis is based on the use of self-similar variables, suggested by 
the invariance of the Burgers equation under the group of scaling transformations 
{x,t,u) I— i- {cx,c^t,u/c) (for subsequent contributions in the same direction, see 
[12]). More recently, it has been shown in [3] that the slow motion for the viscous 
Burgers equation on the whole real line is determined by the presence of a one- 
dimensional center manifold of steady states for the equation in the self-similar vari- 
able (corresponding to the diffusion waves) and a relative family of one-dimensional 
global attractive invariant manifold. In a short-time scale, the solution approaches 
one of the attractive manifolds and remains close to it in a long-time scale. 

At the present day, results relative to metastability in the case of systems ap- 
pear to be rare. Slow dynamics analysis for systems of conservation laws have been 
considered in [11] , basic model examples being the Navier-Stokes equations of com- 
pressible viscous heat conductive fluid and the Keyfitz-Kranzer system, arising in 
elasticity. The approach is based on asymptotic expansions and consists in deriving 
appropriate limiting equations for the leading order terms, in the case of periodic 
data. In [15], the problem of proving convergence to a stationary solution for a 
system of conservation laws with viscosity is addressed, with an approach based on 
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a detailed analysis of the linearized operator at the steady state. A recent contribu- 
tion is the reference [3], where the authors consider the Saint- Venant equations for 
shallow water and, precisely, the phenomenon of formation of roll-waves. The ap- 
proach merges together analytical techniques and numerical results to present some 
intriguing properties relative to the dynamics of solitary wave pulses. 

Summing up, apart for the formal expansions methods, the rigorous approaches 
used in the literature are largely based on typical scalar equations features. The 
first of these properties is the direct link between the scalar Burgers equation and 
the heat equation given by the Hopf-Cole transformation: u = —2e (f)~^dx(l), and the 



consequent invariance of equation (1.1) under the group of scaling transformations 
{x,t,u) I-)- {cx,c^t,u/c). On the one hand, the presence of such a connection is 
an evident advantage, since it permits to determine optimal descriptions for the 
behavior under study (see [131 US EH]); on the other hand, to use such exceptional 
property makes the approach very stiff and difficult to apply to more general cases. 
A different "scalar hallmark" is the base of the approach considered in [7j, where 
the authors make wide use of maximum principle and comparison arguments, taking 
benefit from the fact that the equation is second-order parabolic. 

In order to extend the results to more general settings and specifically for systems 
of PDEs, it is useful to determine strategies and techniques that are more fiexible, 
paying, if necessary, the price of a less accurate description of the dynamics. A 
contribution in this direction has been given in [23], where the location of the shock 
transition for a scalar conservation law in the quarter plane has been proved by 
means of weighted energy estimates, extending the result proved in [19], that used 
an explicit formula -determined by means of the Hopf-Cole transformation- for the 
Green function of the linearization at the shock profile of the Burgers equation. 

The essence of this article is to contribute to the definition of a fiexible language, 
hopefully relevant to more general contexts and, mainly, in the case of systems. 
With this direction in mind, we follow an approach that it is strictly related with 
the projection method considered in [HI EZ] and we go behind the philosophy tracked 
in the analysis of stability of viscous shock waves by K.Zumbrun and co-authors 
(see [nH [221 12I])- Precisely, once a set of reference states {U'^{-;^)} is chosen, we 
separate two distinct phases: 

- to determine spectral properties of the linearized operator at such states; 

- to show that appropriate assumptions on the structure of the eigenvalues of such 
operators together with a control on how far is the approximate state from being an 
exact solution imply the presence of a metastable behavior. 

The main advantage in such a separation stems from the fact that, in principle, 
it should be possible to obtain numerical evidence of special spectral structures in 
cases where analytical results appear to be not achievable. 
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With respect to the framework of shock waves stabihty analysis, there are two 
main differences. First of all, we concentrate on the case of bounded domains and, 
therefore, the spectrum of the linearized operator is discrete and given by a diver- 
gent sequence of (real) values. Additionally, the reference states If^ generically are 
approximate solutions, in the sense that they satisfy the steady state equation with 
an error, that converges to zero as e ^ 0^. Hence the perturbations of such states 
satisfy at first order a non-homogeneous linear equation, with forcing term, formally 
negligible as e — )■ 0"*". 

Our approach consists in approximating the evolution of the couple where 
^ denotes the parameter for the approximate manifold and v the perturbation of 
the profile If^, by a partial linearization, giving raise to a system which we call 
quasi-linear system. This is obtained by linearizing with respect to v and keeping 
the nonlinear dependence on ^, in order to keep track of the nonlinear evolution 
along the approximate manifold. The main contribution of the paper is Theorem 



3.1[ stating that, assuming a number of assumptions relative to the elements of the 
approximate manifold and the linearized operator at such states, and a coupling 
condition, linking the first eigenvalue of the linearized operator with the nonlinear 
operator evaluated at U^, the solution to the quasi-linearized system is described by 
the evolution of a reduced system where the equation for ^ is decoupled from w and 
hence solvable by means of the standard separation of variable method. Details on 
assumptions and statement are given in Section [3j 

We close the Introduction with an overview on the structure of the article. In 
Section [2| we introduce a general framework to describe the dynamics in a cilyn- 
drical neighborhood of an approximate manifold {U'^{-;C,)} of solutions for scalar 
evolution equations, using as coordinates the parameter ^, describing the manifold, 
and a distance vector v, determined by the difference between the solution u and 
the element of the manifold. The couple (^, v) turns to solve an ODE-PDE cou- 
pled system of equations and, in Section |3j we concentrate on the analysis of an 
approximation of the system, obtained by disregaring o{v) — terms, which we call 
quasi-linearized system, since it still preserves nonlinear dependence on ^. In Section 
|4| we analyze spectral properties of the diffusion-transport linear operator, arising 
from the linearization at the state f/^(-;^). Specifically, we introduce our basic as- 
sumptions needed to state and prove Theorem We show that, under appropriate 
assumption on the limiting behavior of as e — ?■ O"*", the spectrum can be decom- 
posed into two parts: the first eigenvalue turns of order ©(e"*"/^), C > 0, hence 
small as £ — )■ 0"*"; all of the remaining eigenvalues are less than —C/e, C > 0. Such 
estimates will translate in the following Sections into a one-dimensional dynamics, 
since all of the components relative all of the eigenvector except the first will have 
a very fast decay for e small, and in a slow motion along the approximate manifold, 
as a consequence of the size estimate for the first eigenvalue. Finally, in Appendix 
(Section [sj, we prove that in the purely hyperbolic case, i.e. £ = 0, the entropic 
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solution to the initial-boundary value problem converge in finite time to one of the 
jump configurations f/^ypl'^O some ^. 

2. General framework 
Let us consider the Cauchy problem for a general evolution equation 
(2.1) dtu = j^''[u], u\^^^ = uo 

where J-"^ denote a nonlinear differential operator, depending singularly on the pa- 
rameter e > 0, so that the formal limiting problem dfU = J^^[u] is of lower order. 



Typically, equation (2.1) is complemented with appropriate boundary conditions. 



appearing in the definition of an appropriate Hilbert space X such that a solution 



to (2.1 ) is a function u : [0, +oo) — )■ X. 



Denoting by m"^ = u^(x, t) the solution of (2.1 ), we are interested in describing the 
behavior of for small e. 

The basic example we have in mind is the initial-boundary value problem for 
a scalar conservation law with viscosity with Dirichlet boundary condition in the 
bounded interval I = {—£,£), that is 

dtu + d^f{u) =edlu xel, t>0 
(2.2) { u{±£,t)=u± t>0 

u{x, 0) = Uq{x) X G I 



for some e,i > 0, u± G M and flux function /, chosen so that assumptions (1.5) 
hold. 

We assume to have a one-parameter family {U'^{-;^)} in X, parametrized by 
^ G /, whose elements are approximate stationary solution to the problem, i.e. 
J^^[f/^(-; ^)] — )■ as £ — i- 0. Precisely, we assume that the term belongs to 

the dual space of the continuous functions space C{I) and there exists a family of 
smooth positive functions Q'^ = ^^{C,), uniformly convergent to zero as £ — )■ 0, such 
that, for any ^ G /, there holds 

m-),j^^[u%;0])\ < I^Lo. G c{i). 

The family {U'^{x; Ol^e/ will be referred to as an approximate invariant manifold with 



respect to the flow determined by (2.1) in the Hilbert space X. 

The dependence of on e plays a fundamental role, since it drives the departure 
from the approximate invariant manifold. In the specific case under consideration, 
such term is exponentially small, meaning that it behaves as e"*"/^, with C > 0. 

Let us stress that, differently to the construction in [16j and in ^27J, where the 
approximate solutions satisfies exactly the equation and the boundary condition to 
within exponentially small terms, here we assume the elements f/^ to satisfy the 
boundary conditions exactly and the equation approximately. 

Example 2.1. In the case of Burgers equation, i.e. /(s) = | and -u^ := =Fw*, 
for some > 0, we consider a function obtained by matching two different steady 
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states satisfying, respectively, the left and the right boundary condition together 
with the request f/(0 = 0; i^i formulas, 

/t_ tanh — x)/2£:) in 

tanh — x)/2£:) in 

where k± are chosen so that the boundary conditions are satisfied 

(2.3) /s:±tanh(|^(^T£)) = n±. 

By direct substitution, denoting by 6x=£^ the usual Dirac's delta distribution concen- 
trated at X = ^, we obtain the identity 

in the sense of distributions. In particular, ?7'^(-, ^) is a stationary solution if and only 
if ^ = 0, that corresponds to the equilibrium position for the shock layer location in 
the case of the Burgers flux. Going further, by differentiation, we have 

In order to determine the behavior of fi^ for small e, we need an asymptotic descrip- 
tion of the values k±. To this aim, let us set k± := =F^^±(1 + h±), so that, denoting 
by A-t := £ =F ^ the distance from ^ to ±i, relation ( 2.3[ ) becomes 

tanhf^^a + Z^i)^ ^ 



Therefore, the values h± are both positive and thus 

tanh ( ^^^±^1 < ^ 



that gives the asymptotic representation 

h± < TT — \ ^ - 1 = A 7 = 2 e±"±^±/" + l.o.t., 

- tanh(T'u±A±/2£) eT"±^±/^ - 1 

where l.o.t. denotes lower order terms. Thus, we infer 



e 

showing that the term [c^^f^'^L^j is null at ^ = and exponentially small for e — )■ 0^. 

Example 2.2. A different example of approximated manifold for a Burgers type 
equation can be obtained by following the traveling wave approach considered in [7] . 
In this case, we look for a solution to (2.2) of the form 

(2.4) u{x,t) = $"(x - V'iOt), satisfying <l>"(±£) = ^m,, $"(0 = 



10 C.MASCIA AND M.STRANI 

for some ^ G {—£,£) and for some u^, > 0. Such solution moves with speed V and 
satisfies the boundary condition only if t = 0. Then the function $^ has to satisfy 

(2.5) edl<^' + {<^' + V)d^<l>' = 0, <^%±i) = ^u. 



Since = (9i$^ by Q we get 

In order to determine the behavior of we need an asymptotic description of 

the speed V in terms of e. As proved in for any ^ G {—i, i) and any e > there 
exists a unique solution of (2.5), satisfying for e — )■ the asymptotic 

^%x) = Tu* + O {exp{-u^\^ - x\/e) + Re) for x^^ 

V%0 = «4exp {-u*{i - 0/e) + exp {-u,{i + 0/^)] + 0{Rl/e) 

where Rs is defined by 

R, := exp {-u^{£ - + exp {-u^i + 0/^) 

Thus we have an asymptotic representation for the speed V 

V ~ M,e-"*(^+^)/" - M,e~"*(^-«)/^ + l.o.t. 

which coincides, up to a multiplicative factor, with the formula determined in Ex- 



ample 2.1 



Once the one-parameter family {U'^{-;C,)} is chosen, we write the solution to the 
initial value problem (2.1) as 

with ^ = C,{t) G / and v = v{-,t) G X to be determined. Substituting, we obtain 

(2.6) dtv = Clv + T^[U%-,0]-d^U%-,0 § + Q^v,^] 

where 

Q^[v, e] := (■; 0+v]- J^nU'i-, 0] - dJ^^lU^-; 0] v. 
Next, let us assume that, for any ^, the linear operator £| has an increasing sequence 
of real eigenvalues A| = A|.(^) with Afc— )■— ooasA;— J'+cxd with corresponding right 
eigenfunctions (pl = 0|.(-;O- Denoting by ipl = 'i/'KsO eigenfunctions of the 
corresponding adjoint operator and setting 

we can use the degree of freedom we still have in the choice of the couple {v, ^) in 
such a way that component f i is identically zero, that is 

I {rA-, m),v{; t)) = and ^o),Vo{-))) = 0. 
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Using equation (2.6), we infer 

+ (■; 0] - d^U^i-; of + Q1^, e]) + {d^^, of , v) = 

Since {ipl, C^v) = Ai(^/'f,f), we obtain a scalar differential equation for the variable 
^, describing the reduced dynamics along the approximate manifold, that is 



(2.7) 
where 



am--=m-,0,d^U%-,0) and a%^,v) := a^O - {d^i'ti-,0,v), 
together with the condition on the initial datum 



To rewrite equation (2.7) in normal form in the regime of small v, we assume 

l«o(OI = \m-;0,d^U'i-.m > Co > G /. 

for some Cq > 0. Such assumption gives a (weak) restriction on the choice of the 
members of the family {U^} asking for the manifold to be never transversal to the 
first eigenfunction of the corresponding linearized operator. From now on, thanks 
to the previous hypothesis, we can renormalize the eigefunction ipf so that 

Since we are interested in the regime f — t- 0, we expand as 



1 + 



+ o{\v\) = l + {d^ijl,v)+o{\v\). 



Inserting in (2.7), we end up with the nonlinear equation for ^ 



(2.8) 

where 
(2.9) 

(2.10) 



dt 



r(0(l + (95^f,t;))+Pie,' 



Using (2.8), the equation (|2.6|) for v can be rephrased as 
(2.11) 
where 



Qnv,^]-d^u^i-,0pn^,v]. 
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Let US stress that, by definition, there holds 

(2.12) {i;ti-,O,H'i-,O) = 0, 

so that H'^{-; ^) is the projection of J-'^[f/^(-; ^)] onto the space orthogonal to 0f (■; ^). 

To show how such formulas can be handled, at least formally, in concrete situa- 
tions, let us analyze problem (2.2), namely 

T'[u]=edlu-dJ{u) 



for e > and / satisfying assumptions (1.5). Retracing the definitions introduced 
above and setting a^(x;^) := /'{U'^^x ;(,)), we get the following expressions 

where the adjoint operator has to be considered with Dirichlet boundary con- 
ditions, and 

where A^^ = o(|f |), so that = o{\v\, \dxv\). Formally, for small e and small v, the 
dynamics of the parameter ^ is approximately given by 

with 6'^ given in (2.9). Next, we need to identify the functions ipl and d^lf^ in the 
limiting regime e — )■ 0, at least approximately. For e ~ 0, the function ipf is close to 
the eigenfunction of the operator relative to the eigenvalue A = 0, with 

Hence, we obtain the representation formula 

1 - e"+(^-«)/^)(l - e-«-(^+^')/=) X < 

(1 - e-"-(^+«)/")(l - e"+(^-^')/=) X > 

so that tpl 1, provided ^ is bounded away from the boundaries ±i. Additionally, 
with the approximation U^{x\^) ~ [/^^^(a;;^), defined in (1.6), we infer 

U^{x-i + h)-U-{x-0 1 
1 

so that we expect d^U^ to converge to — \u\ 5^ as e — )■ in the sense of distributions, 
so that ao(0 ~ ^[""l- Therefore, we deduce an (approximate) expression for the 
function 9^ 

that, with the choice of proposed in Example |2.1[ reduces to 

(2.13) r (0 ~ ^ (e-"*(^+«)A _ 

which coincides with the corresponding formula determined in |27j . 
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3. Analysis of the quasi-linearized system 



Next, let us go back to the system (2.8)-(2.11) for the couple and let us 

neglect the o{v) order terms: 

(3.1) I ft=(^'(0{i + {dcr,M), 

to be complemented with initial conditions 

(3.2) C(0) = Co e i-i, i) and w{x, 0) = wo{x) G L^(/). 



From now on, we will refer to this system as the quasi-linearization of (2.8 )-(2.11 ). 
We are interested in describing the behavior of the solution to such system in the 
regime of small e. 

Shortly, the quasi-linearized system is determined by an appropriate combination 
of the term J-'^[?7^], measuring how far is the function f/^ from being a stationary 



solution, and the linear operator £|, controlling at first order how solutions to (2.1) 
depart from If^ when the latter is taken as initial datum. To state our first result, 
we need to precise the assumption on such terms. 

HI. The family {f/^(-,0} is such that belongs to the dual space of C(J) 

and there exists a family of smooth positive functions fl^ such that 

i(v^(-),-F^r(-,o])i < ^^(0 i^L. G c(/). 

with converging to zero as e — )■ 0, uniformly with respect to ^. 

H2. Let {■ ■ • < A|(^) < • ■ ■ < A|(^) < Af (^)} be the sequence of eigenvalues of 
the linear operator £|. Assume that for any ^ G {—£,£) there hold 

m)-m)>c, m)<o xm<-ck^ withfc>2. 

for some constant C > independent on k, e and ^. 

H3. Given ^ G /, let 0fc(-;O and V^KsO be a sequence of eigenfunction for the 
^ and 



operators C\ and respectively, normalized so that 



(3.3) {rA-, 0, d^U^-- 0) = 1 and (^|, 0|) 



1 ifj = A:, 

iij^K 



Then we assume 

(3.4) = Y^^ruA^^;)' < c wk. 

for some constant C independent on the parameter ^. 
For later use, note that, by differentiation, there holds 

(3.5) {d^rj,<l>l) + {rj,d^<Pl) = 0. 
Also, we use the notation A| := supA|(^). 
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Theorem 3.1. Let hypotheses Hl-2-3 be satisfied. Additionally, assume that 

(3.6) n%0<C\Xim €(-£,£) 

for some constant C > independent on e and ^. 

Then, denoted by {C,w) the solution to the initial-value problem (3.1) -(3.2), for 
any e sufficiently small, there exists a time such that for any t <T^ the solution 
w can be represented as 

w = z + R 

where z is defined by 

z{x,t) ■.= J2wk{0)exp (^J^ Xl{C{a))da^ M^'Xit)), 

and the remainder R satisfies the estimate 

(3.7) \R\^, <Cm,^ jexp Q\l{aa))da^ \wo\^, + l| 

for some constant C > 0. 

Moreover, for wo sufficiently small in L"^, the final time can be chosen of the 
order —C \ In I fi^ Loo- 



The conclusion of the proof of Theorem 3.1 is based on the following version of a 
standard nonlinear iteration argument. 

Lemma 3.2. Let f = f{t),g = g{t) and h = h{s,t) be continuous functions for 
t G [0, T] for some T > 0, such that 

fit) > 0, g(t) > 0, g decreasing, h{s,t) > 0. 

Let y = y{t) be a non-negative function satisfying the estimate 



y{t)< / {f{s)g{t)y\s) + h{s,t)} ds 
Jo 

for any t <T. Then, if, for any t G [0,T] there holds 

(3.8) sup / g'^{s)f{s)ds ■ sup g^^{t) / h{s,t)ds<- 

te[o,T]Jo te[o,T] Jo 4 

then 

y{t) < 2 sup / h{s,T)ds 
Te[o,t] Jo 

for any t G [0, T]. 



Proof of Lemma \372 . The auxiliary function w{t) := g ^(t) y{t) enjoyes the estimate 

w{t) < [ {a{s)w^{s) + /3{s,t)} ds 
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where a{t) := f{t) g^{t) and l3{s,t) = g~\t)h{s,t). Set 

N{t) := sup w{t). 

re[0,i] 

Then, for any t G [0,T], there holds 

w{t)< (^j^ a{s)d^ N\T) + P{s,t)ds 
and, as a consequence, also 



N{T) < AN^{T) + B 



where 



A = A{T) := sup / a{s) ds, 

tG[0,T] Jo 

Since iV(0) = 0, if 1 - iAB > 0, then 

1 - VI - 4Afi 



B = B{T) := sup / I3{s,t)ds. 

tG[0,T] Jo 



N < 



2B 



2A 



1 + VI - 4:AB 



< 2B. 



In term of y, if (3.8) holds, then 



y{t)<2g{t) sup (^-^(t) / h{s,t)ds. 

ie[0,T] Jo 



The final estimate follows from the monotonicity of the function g. 



Lemma 3.2 gives the final step needed to prove the Theorem. 



Proof of Theorem \3.1\ Setting 

j 

we obtain an infinite-dimensional differential system for the coefficients Wj 

dwk 



(3.9) 



dt 



where, omitting the dependencies for shortness, 

F:=H^ + Y, {mi <P] - d^ct>) §} = - E («^- + E 

3 3 I 

and the coefficients a^, are given by 

Convergence of the series is guaranteed by assumption (3.4). 
By (3.5), for the coefficients there hold 
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SO that we can also take advantage from the relation (^/'f , aj) =0 for any j. Thanks 
to these relations, equation (3.9) for k = 1 simplifies to 

(3.10) ^ = XliO w, - r (C) b,e) w, wj 



Now let us set 



Note that, for < s < t, there hold 

Ek{s,t) = ^!^ and 0<E,{s,t)<e^^^'-'\ 
Ek{0,s) 

From equalities ( 3.10[ ) and ( |3.9[ ), choosing wi{0) = 0, there follow 
wi{t) = - [ e%C)Y,{^lbje)wiWjEi{s,t)ds 
Wk{t) = Wk{0) Ek{0,t) 

for k > 2. Such expressions suggest to introduce the function 

z{x,t) := J2M0)Ek{0,t)<Pl{x;C{t)), 

k>2 



which satisfies the estimate \z\^2 ^ l^o|^2^ 

From the representation formulas for the coefficients Wk, since 

\9%C)\<cn%C) and m,H^)\<cn%0{i + m,d^u^)\} 

for some constant C > depending on the norm of ipl, there holds 

2 



\w-z\l^<c[f Q%0 J2m, 5^05)1 HI 5^ K^^V^L 0^)1 \we\ E,{s,t) ds 

k>2 ^-^0 ^ 

j j e 

<cU'Q%0\M%Ei{s,t)dsy + Cj2{f^%0{^ + \M%)E,{s,t) ds 

Jo 7,^0 ^0 



METASTABILITY FOR SCALAR CONSERVATION LAWS 

Since \/a + b < ^J~a + we infer 

<C fn%0\w\\E,{s,t)ds + Cy fn%0{l + \w\\)Ek{s,t)ds 
Jo ^ t^io 

<C f n'iO{\w\l,E,is,t) + {l + \w\'J J2Ekis,t)}ds. 

"^0 fc>2 
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\w — z\ 



The assumption on the asymptotic behavior of the eigenvalues can now be used 
to bound the series. Indeed, there holds 

E,{s, t) < E,is, t) Y: <C{t- E,{s, t) 

k>2 k>2 ' ^ 

As a consequence, for unknown w such that {wl^^ — ^ some M > 0, we infer 
Ei(t,0)|w - z\^^ < C j^n%0[\w - z\l^ E,{s,0) 

+ \z\l^ Ei{s, 0) + (t - ^^^^^ 0)} ds 



Let us set 



N{t) := sup \w — z\^^ Ei(s,0) 



Then, since A| < Af , we obtain 
Ei{t,0)\w - z\^^ <C n%ON\s)Ei{0,s)ds 



+ 



By assumption ( |3.6 ), Af < — Cl)'^ for some constant C > 0, hence 



n%()N\s)Ei{0,s)ds< / n%()N\s) exp i -C / Q^Qda] ds 



Moreover, there holds 



g(2A|-Af ). < / e^'^'ds 







[e^2« _ 1) < 



|A 



[\t - s)-^/^ E2{s, t) ds < fit - sY^I^ e 
Jo Jo 

so that, recalling that A| is bounded away from 0, 
Ei(t,0)|w - < c|iV=^(t) l-exp^^-C j Q'iOda 



|A||V2 



+ C\Q^L {\wo\l,+E,{t,0))] 
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Next, we end up with the estimate 



N{t) < AN'^{t) + B with 



Hence, as soon as 



A:=c|l-exp (-C^ ^'iOda 
B ■.= C\n'\^^(\wo\l^+E,{t,0) 



(3.11) AAB = AC^\n'\^^(^\wo\%+E,{t,0)'^(^l-exp(^-C J^n^iOda^^ < 1 



there holds 



N{t)< 



2B 



<2B = C\n'\^^{\wo\l^+E,{t,0) 



1 + Vaab 

that means, in term of the difference w — z, 

\w-z\^, <C\n'\^^(\wo\l^E,{0,t) + l^ 

Condition (3.11) gives a constraint on the final time T^. Since 1 — e~^-fo ^HO^a- ^ 
it is enough to ask 

(3.12) AC^ (^\wo\l^ + E,{t, 0)) < 1 



to assure condition (3.11) is satisfied. Constraint (3.12) can be rewritten as 

C 



C exp ^]"(C) < exp da^ = Ei{t,0) < 



2 

l2' 



SO that we can choose of the form 



In 



C 



for wq sufficiently small. 



□ 



As a consequence of the estimate (3.7), for \w\^^ < M for some M > 0, the 



function ( satisfies 



dt 



^^(C)(l + r) with \r\<C{\wo\^,e'''^' + \n\^). 



where the constant C depends also on M . In particular, if e and \ wq\ ^ are sufficiently 
small, the function C, = C,{t) has similar decay properties with respect to the function 
7], solution to the reduced Cauchy problem 

dt] 



dt 



This preludes to the following consequence of Theorem 3.1 
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Corollary 3.3. Let hypotheses Hl-2-3 and (3.6) be satisfied. Assume also 

r (0 < 0. 



(3.13) 



s 6'^{s) < for any s G /, s 7^ 



and 



Then, for e and |ifo|^2 sufficiently small, the estimate (3.7) holds globally in time 
and the solution (C,"?^) converges to (CO) as t ^ +00. 



Proof. Thanks to assumption HI, for e and |wo|^2 sufficiently small, estimate (3.7) 
holds. Hence, for any initial datum Co, the variable C, = C{t) solves an equation of 
the form 

^ = r(C)(l + r(t)) with |r(t)|<C(KL2e^^* 



Therefore, ({t) converges to ( as t ^ +00 and the convergence is exponential, in 
the sense that there exists < such that 



(3.14) 



|C(t)-Cl<ICo|e^^ /3^~r(C) 



for any t under consideration. Furthermore, from (3.9), we deduce 



Wk{t) = Wk{0) exp Xlda j + J {^IjI,F){s) exp ( / Xlda ] ds 
Thus, by the Jensen's inequality, we infer the estimate 



w\%{t)<C\\wo 



12 2Aft 



l2 



[s e 



2Af(t-s) 



Let z/^ > be such that (t) < Ce"'' *; then, if 7^ |Af|, there holds 

\w\l^ {t)<c{ \wo\% e^^t * + t [e-^^' * + e^^^ *) } 
showing the exponential convergence to of the component w. 



□ 



Estimate (3.14) shows the exponentially slow motion of the shock layer for small e. 
Precisely, the evolution of the location of the shock towards the equilibrium position 
is much slower as e becomes smaller, since — ?■ as e — ?■ 0. For example, when 
/(s) = sV2, C = and ^'(0) ~ 6"^/^ (see formula ( |2l3| ). 

Let us also stress that in the regime (C, vo) ~ (C, 0), a linearization at the equilib- 
rium solution f/^(x;C) would furnish a more detailed description of the dynamics, 
since the source term due to the approximation at an approximate steady state 
would not be present. In fact, the description given by the quasi- linearization is 
meaningful in the regime far from equilibrium and its aim is to describe the slow 
motion around a manifold of approximate solutions. 
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4. Spectral analysis for the diffusion-transport operator 

Our concern in the present Section is to establish a precise description on the 
location of the eigenvalues of the linearized operator, in order to show that the 
general procedure developed in the previous Sections is indeed applicable in the 
case of Burgers equation. The problem of determining the limiting structure of the 
spectrum of the type of second order differential operators we deal with has been 
widely considered in the literature. Among others, let us quote the approach, based 
on the use of Priifer transform, used in [5] , in the context of metastability analysis 
for the AUen-Cahn equation. Here, we prefer to follow the strategy implemented 
in [13], for the linearization at the steady state of the Burgers equation. In what 
follows, we show that the same kind of eigenvalues distribution holds in a much 
more general situation, the main ingredient being the resemblance of the coefficient 
a"^ to a step function a", jumping from a positive to a negative value, as £ — )■ 0"*". 



Fixed e > and linearizing the scalar conservation law (2.2) at a given a reference 
profile f/^ = [/'^(x;^), satisfying the boundary conditions f/^(±£;^) = u±, we end 
up with the differential linear diffusion-transport operator 

(4.1) C^u := e dlu - d^ia' u) u{±i)=0, 

where = a^(x;^) := /'(f/^(x; ^)). The aim of this Section is to describe the 
structure of the spectrum o"(£|) of the operator £^ for e sufficiently small. 
Given the function a'^, let us introduce the self-adjoint operator 

(4.2) M^v := d^v -Wv v{±e) = 0, 
where 

1 \^ 1 da^ 



(^■^' ''■■=[2"') ^^d.- 

By omitting the dependencies from ^ for shortness, a straightforward calculation 



shows that if u is an eigenfunction of (4.1) relative to the eigenvalue A, then the 
function v{x) defined by 

v{x) = exp (^-^ j (^'(y) dy^ «W 

(with xo arbitrarily chosen) is an eigenfunction of the operator A^^ relative to the 
eigenvalue /i := eX. Since is self-adjoint, we can state that the spectrum of the 
operator £^ is composed by real eigenvalues. Moreover, if u is an eigenfunction of 
(4.1) relative to the first eigenvalue Af, integrating in the relation C^^u = Am, 



we deduce the identity 

0=1 {C\-\\)udx = e{d^u{t)-d^u{-l)) -\\ I u{x)dx 
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Assuming, without loss of generality, u to be strictly positive in and normal- 

ized so that its integral in is equal to 1, we get 

XI = e {d,u{l) - d,u{-£)) < 

Hence, for any choice of the function a^, there holds 

a(£|)c(-oo,0). 

Our next aim is to show that under appropriate assumption on the behavior of the 
family of functions a*^ as £ ^ 0"*", it is possible to furnish a detailed representation of 
the eigenvalue distributions for small e. Specifically, we are interested in coefficients 
behaving, in the limit £ — > 0"^ as a step function of the form 

a_ xe(-^,0, 



(4.4) a%x) := 



for some ^ G (— ^, ^) and a+ < < a_. We will show that, under appropriate 
assumptions making precise in which sense "resemble" for e small, the first 
eigenvalue Af turns to be "very close" to for e small, and all of the others eige- 
navalues A|, with k >2, are such that £A| = 0(1) as £ — > O'*'. 

Estimate from below for the first eigenvalue. We estimate the first eigenvalue 
III of the operator Al| by means of the inequality 



1^1, 



for smooth test function ip such that = 0. Let us consider as test function 

tp~{x) := iI'q{x) — K^{x), where 

V'o(^) — exp ^ a'{y) dy^ , 

A direct calculation shows that M^^t/^ :— If K and, assuming the family to be 
uniformly bounded, we infer 

\h^K'\ , \K'\ , C 





\K^\ 


i2 




L-2 ~ 


\K^\ 


1,2 



lz,2- 



as soon as \'4'l\^2 > 

The opposite case being similar, let us assume ■0o(— ^) > V'o(^)- Prom the defini- 
tion of K^, it follows 

I^'IL = ^{V^o(^) + V^o(^)^o(-^) +^o'(-^)} < S^V'o'l-^)- 
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Therefore, we deduce 



where 



Since converges to the step function as e ^ 0^, it is natural to approximate 
the latter integral in term of the corresponding one for a°: 

r = j' exp j\a' - a!'){y) dy^ exp Q ^ a%y) dy^ dx > e'^^^-^^^L^^' jO. 



Since, for e small. 



^a-{x+£)/e _|_ ga_(^+£)/£ / ^a+{x-^)/e 



I a_ a+ J a_a+ 

the subsequent estimate holds 

whenever — a°|^^ < cqE for some cq > 0. Thus, we deduce for the first eigenvalue 
of the self-adjoint operator the estimate | < Ci e^'^^^ for some positive con- 
stant Ci, C2. As a consequence, since the spectrum o"(£|) coincides with £~^(j(A^|), 
the next result holds. 

Proposition 4.1. Let be a family of functions satisfying the assumption: 
AO. there exists C > 0, indipendent on e > 0, such that 

da^ 



\^ 1 1^00 ~t~ s 



dx 



< C 



If there exists ^ G a+ < < a_ and C > /or which — a^l^i < Ce, t/ien 

there exist constants C, c > swc/i that —C e^^^^ < < 0. 

Let us stress that the request a+ < < a_ is essential, even if hided in the proof. 
If this is not the case, the term would not be small as e — )• O"*" and its norm 
would not be bounded by the L^-norm of ip^. In fact, the statement in Proposition 



4.1 may not hold when a± have the same sign, the easiest example being the case 
a"^ = a+ = a_ > 0. 

The next Example gives a heuristic estimate for the first eigenvalue 

Example 4.2. Given —a < < /3 and a± G M, let us set / = (—a,/?), [a] := a+ — a_ 
and 

a{x) = a_X(_,_o)(a;) + a+x^o.^^{x). 
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Given A"^ > 0, let us look for function u E C{I), such that 

Cu := eu" + (a(x) u)' + \^ u = 0, u{—a) = u{(3) = 

in the sense of distributions. Since a' = [a] 6q, this amounts in finding two functions 
such that 

C±u := eu'^ + a± + u = 0, a) = = 

and the following transmission conditions are satisfied 

M+(0) - M_(0) = and e (m'+(0) - m'_(0)) + [a] m±(0) = 0. 
The characteristic polinomial of C± is A"^) := e ^"^ + a± ^ + with roots 
— a_ ± A_ I —a^ ± A_i_ 



± 



2e 



2e 



where A± := ^/(a±y^^^^F}^ . 



Assume A^ < {a±y/Ae. Choosing u± in the form 

M_ {x) = A_ {e"- ("+^) - e^- ("+^) ) and m+ (x) = (e'^^ ^^-"^^ - e'^^ ('^"^^ ) . 
Setting := e^-" and := e'^'^^, there holds 

u_{o) = A^{0+ - r ) u_{o) = A^ifitel - iizez) 

w+(0) = A+{eX - 9-) <(0) = - fi-9-). 

Therefore, the transmission conditions take the form of a linear system in A± 
' {9% - 9-^)A^ - {9t - r )A_ = 0, 

(4.5) ] {(2^/"t + H) K - i^e^i- + [a]) 9-}a+ 

+ I {-26 /i+ + [a]) 9t + {2e fiZ - [a]) 9z^A_ = 0. 

After some manipulation, the determinant D = D{X'^,e) of (4.5) can be rewritten 
as 

D = {[a] + [A]) 9X91 - (H - {A}) 9-9t - ([a] + {A}) 9X9Z + {[a] - [A]) , 

where [A] := A+ - A_ and {A} := A+ + A_. 

Since — 4 X = \k\ — 2\k\^^ x + o(x), for — )■ 0, there hold 



e\n{9X9^_ 



\a_\a 



/3 



e ln(6'_^6'l) = (a+ /3 + |a_ |a) - 



/3 



£:A^ + o(eA 



+ 



eA + o(£A), 



£ln(^+r 

eln(^;r 



a_|_ 
a+/3 - 



a 

a_ I 

A 



eX' + o(£A") 



a 
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{A} = ^{a+Y -Ae\' + ^{a.f - Ae\' = [a] 1 + + o{e\ 



[A] = ^/{a+y-Ae\' - ^/{a^y-Ae\' = {«} ( 1 - 
Hence, we infer 



_g(a+/3+|a_|a)/£ _ Q_ga+^/e 



so that D ~ for 



£ a 



For a± = ±u^, a = i + ^ and (3 = i — ^, the above expression becomes 



A- ^ ;i (e-*(^-OA + e-«.(^+C)A) = !^ cosh(n4/£) e" 



Ye ^ ■ e 

to be compared with the expression for Vf^ obtained in Example |2.1 



2m? 



Aul 



sinh(M^, i/e) \ e 



SO that 



A^ 



4|tanh(M*^/£:)| < 4 



Let us stress that this formula shows that hypothesis (3.6) is verified heuristically 
for Burgers type equations. 



Estimate from above for the second eigenvalue. Controlling the location of 
the second (and subsequent) eigenvalue needs much more care and, also, a number 
of additional assumption on the limiting behavior of the function as £ ^ 0+. 
Precisely, we suppose G C°([— £, £]) satisfies the following hypotheses: 
Al. the function a*^ is twice differentiable at any x ^ ^ and 



< < a" 



m 



• 0, 



and 



da'^ 



d^a' 



< < 



dx ' dx"^ ' ' ^ dx dx"^ 

A2. for any C > there exists C(, > such that, for any x satisfying |a; — ^| > c^e. 



there holds 



a°| < Ce 



and 



da^ 



dx 



<c- 



A3, there exists the left /right first order derivatives of at ^ and 

da"^ 



liminf e 

£-5-0+ 



dx 



> 



As a consequence, the function + eX"^ satisfies a number of corresponding proper- 
ties, listed in the next statement. 
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Lemma 4.3. Let the family he such that hypotheses Al-2-3 are satisfied, and let 
< &e such that 

1 



(4.6) 



inf eX" > 

£>0 



-a' 
4 » 



where a„ 



min{|a_ I, |a+|}. 



Then there exist > such that, for e < e^, the functions +eX^ , with defined 
in (4.3), enjoy the following properties: 

Bl. the function 6^ + eX^ is decreasing in (— ^) and increasing in (^,^); 

B2. there exist C,c > such that, for any x with |x — ^| > ce there holds 
¥ + eX'>C> 0; 

B3. there exist the left/right limits ofb^ + eX^ at ^ and 

(3 := limsup(6"(^±) + eX') < 0; 

e-s>0+ 

Proof. Property Bl. is an immediate consequence of assumption Al, since 

d 



dx 



1 p da'^ 1 d'^a^ 



From A2, given C > 0, for x < ^ — e, there holds 

b' + eX'> -(a^ + a°)(a"-a°)--£ 
~ 4 2 



da'^ 



dx 



+ eX^ + -a^_ 
4 



>eX^ + \al-\{l + \a'\e+lce' 
- 4 ° 2 



C 



From such inequality, by choosing C > sufficiently small, and combining with an 
analogous estimate on (^ + c£:,£), property B2. follows. 

For what concerns B3, we observe that, since a(^) = and A < 0, there holds 



1 c/a^ 

limsup(6'^(^±) + eA*^) < limsup - e — — (^) 

e^0+ e^0+ 2 dx 



lim inf e 

£-S>0 + 



da^ 
dx 



(e±) 



< 0, 



thanks to A3. 



□ 



For later reference, we denote the zeros of b^ +eX^ , with — £ < < ^ < < 
Since property B2 holds, we deduce that \y^^ — i\ < £. 



Assume the assumption of Lemma 4.3 to hold, and let A| and yu| = e A| be the 
second eigenvalue of the operators £| and A^^, respectively, with corresponding 
eigenfunctions (/)| and V'f- Such eigenfunctions are linked together by the relation 



(4.7) 



V'K^) = ^ 6xp 



1 

Ye 



a'iy) dy (f)l{x 



for some constants A and x^,. Since A| is the second eigenvalue, the functions 0| and 
ip2 possess a single root located at some point Xq G (—£,£). The sign properties of 
+ described in Lemma 4.3 imply that Xq G {yt^y^.)- Then, 0| and ip2 restricted 

are eigenfunctions relative to the first eigenvalue 



to the intervals 



Xn ana x; 



05 ' 
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of the same operator considered in the corresponding intervals and with Dirichlet 
boundary conditions. 

From now on, we drop, for shortness, the dependence on e of X2, (f>2,4'2,Xo, we 
assume, without loss of generality, Xq > ^ and we restrict our attention to the 
interval J = {xQ,i). Integrating on J, we deduce 

A2 / (f)2dx = e (9^02(^) - dcc4>2ixo)) < -ed^(f>2ixo) 

J XQ 



having chosen 02 positive in J. Assuming 1^2 to be given as in (4.7) with A = 1 and 
x^ = Xq, and normalized so that max?/'2 = 1, from the latter inequality we infer the 
inequality 

(4.8) IA2I >£/-'9,V'2(xo), 

where 



/ := y exp j a'{y)dy^ dx 



'XQ JXQ 

Our next aim is to deduce an estimate from above on and an estimate from below 
for dxip2{xo), in order to get a control on the size of the second eigenvalue A2. 
From the definition of Is, since xq > ^, it follows 



whenever — <Ce. Thus, estimate (4.8) provisionally becomes 

(4.9) |A2|>C^(xo) 

for some positive constant C, independent on e. 

Let xm be such that ip2{xM) = 1, minimum with such property. From the as- 
sumption on ¥ + e A, it follows xm & {xq, ?/+). Then there exists x^ E (xq, xm) such 
that 

-1— (a;L) = > T > — . 

ax Xm — Xq y+ — ^ c^e 

Since the function ip is concave in the interval (xq, y+), we deduce 

f2(.o) > > -. 



Plugging into (4.9), we end up with 

(4.10) IA2I > -. 

e 

for some C independent on e. 
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As a consequence, we can state the next result relative to the second eigenvalue 

Proposition 4.4. Let he a family of functions sastisfying Al-2-3 then there exists 
a constant C > such that A| < —C/e for any e sufficiently small. 



Spectral estimates. Collecting the results of Propositions 4J. and AA_ give a 
complete description for the spectrum of operator £| for small 5, under assumptions 
AO- 1-2-3 on the family of functions a^. 

Corollary 4.5. Let he a family of functions satisfying the assumptions AO- 1-2-3 
for some ^ G a+ < < a_. Then there exist constants c,Ci,C2 > such 

that 

K < -Ci/e and - C2 e""/" < A^^ < 0. 

for any k >2. 

Hypotheses AO-1-2-3 are satisfied in the case of a family of function a"^ that is a 
(small) perturbation of a function with the form 

for some decreasing smooth bounded functions A±, bounded together with their 
first and second order derivatives, and such that ^±(±00) = a± and A'_^{±oo) = 0. 

5. Appendix. The hyperbolic dynamics 
In this Section, we concentrate on the dynamics of the scalar conservation law 
(5.1) dtu + dj{u) = 

with X G (— £, i), together with boundary condition u{±i, t) = u±, and a given initial 
datum u{x, 0) = Uq{x). Our aim is to give a self-contained proved of the finite-time 
stahilization of the solution under appropriate assumption on the flux function / 
and on the boundary values u±. This kind of properties has been showed for the 
first time in [IB] in the case of the Cauchy problem. 



Theorem 5.1. Assume the function f to he uniformly convex, i.e. f"{s) > cq > 
for some constant cq. If u^,u+ are such that < < u_ and f{u^) = f{u_), 



then, for any Uq G 5F(— £, £), the solution u to the initial value prohlem (5.1), 

u{±i,t) = u±, u{x,0) = Uo{x) is such that for some T > and ^ G there 
holds 

u{;T) = U,J--0 m (-£,£) 

'^here U^^^{x]^) := M_x,_,,^)(a;) + u+X(j,,)(x). 
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To prove the statement, we use the theory of generalized characteristics, intro- 
duced in [6]. The convexity assumption on the flux function / guarantees that for 
any point {x,t) £ {—£,£) x (0,+oo) there exist a minimal, respectively maximal, 
backward characteristics, which are classical characteristic curves, hence a straight 
lines with slope f'{u{x — 0, t)), resp. f'{u{x + 0, t)). 

The boundary conditions are understood in the sense of Bardos-leRoux-Nedelec 
[2j, meaning that the trace of the solution at the boundary is requested to take 
values in appropriate sets. To be precise, let u^, G be such that f'{u^) = 

and set 

w if 3 w 7^ M s.t. f{w) = f{u), 

Then, skipping the details (see [20]), the boundary conditions u(±£, t) = u± translate 
into 

u{-i + 0,t) G (-oo,7^n_)] U u{i-0,t) G U [7^M+, +oo) 

Since f{u+) = f{uJ), there holds TZu± = u=p, and the condition can be rewritten as 

(5.2) u{-i + 0,t) e {-oo,u+]U {u_}, u{i -0,t) e {u+}U[u_,+oo) 

In particular, characteristic curves entering in the domain from the left side x = —i 
(respectively, from the right side x = i) possess speed f'{uJ) (resp. speed f'{u^)). 



Uu :-- 



Now, we are ready to prove Theorem 5.1 



Proof. Let u = u{x,t) be the solution to the initial-boundary value problem under 
consideration with initial datum uq. For later use, we set 

C-(t):=sup{xGH,£] ■.u{y,t) = l y y e {-i,x)} U -i, 

C+(t) :=inf{xG ■.u{y,t) = -l VyG(x,£)}U£. 

In particular, C_ < C+. We are going to show that C-(^) = C+(^) fo^' some T > 0. 

1. There exists Tq > such that u{x,t) G for any x G {—£,£). 



Indeed, let u = u{x,t) be the solution to the Riemann problem for (5.1) with 
initial datum 

r X < —i, 

I max|M_, sup Moj x > —i, 

Hence, the restriction of u to {—£, £) x (0, oo) is a super-solution to the initial bound- 
ary value problem under consideration and, by comparison principle for entropy so- 
lution, we infer u{x,t) < u{x,t). Since u{x,t) = m_ for any x < f'{uJ)t — i, there 
holds 

u{x, t) < u_ for X G (-£, i),t> 2i/f'{u_). 

A similar estimate from below can be obtained by considering as subsolution the 



restriction of u to (— £, i) x (0, oo), where u is the solution to (5.1 ) with initial datum 

min{M_,., inf Uq} x < i, 

X > i, 



Uo[X 
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Prom now on, we assume that the solution u takes values in the interval [u^, u+]. 

2. Assume that —i < C-{t) < C+{t) < i for any t; then there exists Ti > such 
that u{C,-{t) + 0, < and it+ < u{C,j^{t) — 0, t) for any t > T\. 

If u is continuous at (C-(t),t) for some r > 0, then it(C-(r) + 0,t) = it_. 
Therefore, the maximal backward characteristic from (C_(r),r) is the straight line 
X = C-(t) + f'{'^-){t ~ For T > 2L/ f'{u ), such curve intersects the boundary 
X = —£ at some a G (0, r). By continuity, all of the maximal backward characteris- 
tics from ({, r) with { > C- (t) and suffuciently close to C- (t) intersect the boundary 
X = —£ at some time cr*(^) smaller than a and close to it. Because of the boundary 
conditions, this may happen if and only if u{C,,t) = m_. Hence, u{x,t) = u_ for 
X e (C_(t),C_(t) + e) for some e > 0, in contradiction with the definition of C_. 
Thus, continuity of u at (C-(t),t) may happen only for r < 2L//'(m_). A similar 
assertion holds for C+. 

3. There exist T > and ^ e [-i, i] such that u{x, t) = U^^^{-; ^) for any t > T. 
Given ^ > 0, let Tg := 2£/e be such that 

:= u{C-{Te) + 0, Te) < u_ and u+ < := u{C+{Tg) - 0, Tg). 

Let x^_ be the maximal backward characteristic from ((^-(Tg), Tg), whose equation 
is x = (-{Te) + f'{u^_){t — Tg). If x^_ hits the right boundary x = i at some 
positive time, the solution u coincides with U^^^^{x; (-{Tg)). Otherwise, there holds 
C4Te) - f'{ut)Tg < £, which gives 

^0 J-e 

Similarly, let be the maximal backward characteristic from {C,^{Te) ,Te) , whose 
equation is a; = C+(^e) + f'{u^+)(t — Tg). If does not intersect the left boundary 
X = —i at some positive time, there holds f'{u^^) < 9. 

Hence, for any £ > 0, we can choose sufficiently large so that > — e and 
< + £. Thus, we have 

(K±_<K- ^ f{u+) - f{u^+e) _ f{u-) - f{u^ -e) 
dt dt u-^ — — e u- — + e 

which is uniformly negative for e sufficiently small. Hence, the curves C+ and C_ 
intersect at some finite positive time T > 0. □ 
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